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Abstract 

We examine the stability of the Garfinkle-Horowitz-Strominger (GHS) black hole under charged 
scalar perturbations. We find that different from the neutral scalar field perturbations, only two 
numerical methods, such as the continued fraction method and the asymptotic iteration method, 
can keep high efficiency and accuracy requirements in the frequency domain computations. The 
comparisons of the efficiency between these two methods have also been done. Employing the 
appropriate numerical method, we show that the GHS black hole is always stable against charged 
scalar perturbations. This is different from the result obtained in the de Sitter and Anti-de Sitter 
black holes. Furthermore we argue that in the GHS black hole background there is no amplification 
of the incident charged scalar wave to cause the superradiance, so that the superradiant instability 
cannot exist in this spacetime. 

PACS numbers: 04.50.Kd, 04.70.-s, 04.25.Nx 


1 Introduction 


Perturbation around black holes has been an intriguing subject of discussions in the past three decades. 
This is mainly because that the study of black hole perturbations is a powerful tool to disclose the 
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stability of the black hole spacetime. If the black hole is unstable against small perturbations, it 
will inevitably disappear or transform dynamically into another object. Stability analysis of (3+1)- 
dimensional asymptotically flat black holes, such as Schwarzschild, Reissner-Nordstrom (RN) and 
Kerr black holes, has been studied thoroughly against different kinds of perturbations including neu¬ 
tral scalar, electromagnetic and gravitational perturbations and all these black holes were found stable. 
Considering that our universe may have a small positive cosmological constant, the perturbation anal¬ 
ysis has been extended to (3+l)-dimensional Schwarzschild-de Sitter (dS), RN-dS, and Kerr-dS black 
holes and evidences for the stability of these more possible astrophysical black holes have been re¬ 
ported. Recently, motivated by the discovery of the correspondence between physics in the anti-de 
Sitter (AdS) spacetime and conformal field theory (CFT) on its boundary (AdS/CFT), the perturba¬ 
tions around four-dimensional AdS black holes have been examined and (3+l)-dimensional AdS black 
holes were found stable under neutral scalar, electromagnetic and gravitational perturbations. It was 
concluded that all of the considered four-dimensional black holes tested for stability are stable, except 
the string theory generalization of Kerr-Newman black holes whose stabilities have not been tested 
due to the difficulty in decoupling the angular variables in their perturbation equations. For a review 
on this topic, see mm for example and references therein. 

Recently, the exploration of the black hole spacetime stability has been extended to examine 
the perturbation against the charged scalar field. In the AdS spacetime, it was first observed that 
the (3+l)-dimensional RN-AdS black hole can be destroyed and the four-dimensional AdS hole can 
become unstable due to the condensation of the charged scalar hair onto the black hole. The (3+1)- 
dimensional RN-AdS black hole will finally be transformed into another new hairy black hole under 
the small charged scalar field perturbation 13 II El, see review for example [6]. It. would be interesting 
to ask whether the observed instability only happens for the AdS black holes because of their special 
spacetime properties, whether such dynamical instability can also appear in other four-dimensional 
black hole backgrounds. In [7], a new instability in the four-dimensional RN-dS black holes against 
charged scalar perturbations was disclosed. This result was later confirmed in [8j. Can this instability 
be a general property? In this paper, we would like to examine this problem further. We will 
extend the discussion of the charged scalar perturbation to the (3+l)-dimensional dilaton black hole 
obtained by Garfinkle, Horowitz and Strominger |9](GHS) from the low energy effective action in 
string theory. Its limiting case reduces to the four-dimensional Schwarzschild black hole, see review 
|10j . The stability of this dilaton black hole has been proved by examining the perturbation against 
neutral scalar fields HU E2EE3]. Here we are going to test the stability of such black hole against 
charged scalar perturbation and try to answer whether the instability observed for the charged scalar 
perturbation can also happen in this stringy four-dimensional black hole background and its limiting 
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Schwarzschild spacetime. 

We will concentrate on the frequency domain studies of the charged scalar perturbation around the 
stringy black hole. It is important to calculate the frequency of perturbations with very high accuracy 
because considerable changing of black hole parameters frequently changes perturbation frequency just 
by a few percent. With the nonzero charge of the scalar field, we will explain that not all available 
numerical methods for solving the eigenvalue problem of the perturbation can keep high accuracy. The 
continued fraction method (CFM) [ HI 15] and the asymptotic iteration method (AIM) |16l T7,18] still 
win the accuracy and efficiency competition against the other numerical methods. The CFM method 
was argued as the most accurate in calculating the frequency of the perturbation [T9]. Although this 
still holds when the charge of the scalar field is low, we will show that its accuracy and efficiency will 
be decreased and lower than that of the AIM method when the scalar field is heavily charged. Further 
comparison of the accuracy and efficiency of numerical methods in calculating the frequency of the 
perturbation is important, because this can help to choose the right numerical tool in obtaining the 
real physics on the stability of the black hole. 

The organization of the paper is as follows. In Sec.2 we will introduce the background spacetime 
of the stringy black hole and derive the equation of motion for charged scalar perturbations. In Sec.3, 
we will first review the CFM and AIM methods for numerical computations of the frequency of the 
perturbation. Then we will compare the efficiency and accuracy of these two methods and also with 
other methods. In the following section, we we will give numerical results on the frequency of the 
charged scalar perturbations. We will present our summaries and discussions in the final section. 

2 The GHS black hole and equation of charged scalar perturbation 

The GHS black hole is a solution obtained from the low energy effective action in string theory by 
dropping all the fields except the metric g^, a dilaton 4> and a Maxwell field F llL/ . In string frame, 
the action is |20| 

W= ikJ dAx ^9 e ~ 2<P i R + 4 W) 2 - i>), (1) 

where we have added the perturbing charged scalar field •0 to study its perturbation on the background 
of the GHS black hole. Here D^ipD^ = (<9^ + iqA^)'ijj(d u — iqA u )il> with q being the charge of the 
scalar field if). V(ip,ip) is the potential of the perturbing charged scalar field. Its usual form is taken 
as + j(ip'ip) 2 , in which g is the mass of the scalar field tj] and the A term represents 

the self-interaction of if:. We set A = 0 hereafter for simplicity. This action can be viewed as a special 
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case in scalar-tensor theory [21 j with = e~ 2< ^ and Z((j)) = —4e~ 2< ^. 

Doing a conformal transformation the action can be rewritten in the Einstein frame 

as 

W= lkJ ~ W) 2 - e- 2 ^FH - - e 4< V(</>, ■0). (2) 

At the first sight, it seems difficult to find the exact background solution to equations of motion 
derived from the action Eq. ([2]) (The background solution means the solution not taking into account 
the back-reaction of the perturbing field i/j). However, the symmetry property of this action allows 
one to obtain a one-parameter family of solutions mmm, 


ds 2 E 




+ r 



(3) 


This is the well-known Garfinkle-Horowitz-Strominger (GHS) black hole. Here dll 2 = dO 2 + sin 2 6d(f 2 . 
M is the physical mass and Q the physical charge of the GHS black hole. The electric field and the 
background dilaton field are 
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(4) 


It is obvious that the electric charge and the dilaton are not independent. The GHS solution reduces 
to the Schwarzschild black hole in the limit Q 0. 

In the limit Q —>• 0, the area of the sphere r = Q 2 /M is zero so that this surface is singular. When 
Q 2 < 2AI 2 , this singular surface is surrounded by the event horizon r = 2 M. As one increases Q, the 
singular surface can coincide with the horizon when Q 2 = 2 M 2 and even moves outside the horizon 
and becomes timelike if Q 2 > 2M 2 . In our paper, we will consider the case when Q 2 < 2M 2 . 

Ignoring the self-interaction term, the equation of motion of the perturbing charged scalar field is 



Here g 2 e 2< ^ plays the role of the effective mass square of %). Since the background spacetime is 
spherical symmetry, we can separate the radial and angular part of Using the ansatz if; = 
e~ luli 1 = } = Y(0, ip) and introducing the tortoise coordinate dr = (l — ^-) dr*, we get the radial 

\l r2 -%r r 
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part of the perturbation equation 


in which 
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and the effective potential V e ff = V — q ® . Here l is the spherical harmonic index. It can be shown 
that we always have V > 0 when r > 2 M and Q 2 < 2 M 2 . 


3 Numerical methods 

A practical tool for testing stability of black holes is the numerical investigation of the perturbations 
around black hole backgrounds. Usually considerable changing of black hole parameters results in 
the change of just a few percent in the frequency of the perturbation. Thus the high accuracy of 
the computation is the key factor in examining the perturbation around black holes. Meanwhile 
the efficiency of the computation is also an important factor in solving the perturbation equations 
numerically. 

In this work, we will concentrate on the frequency domain to disclose the property of the per¬ 
turbation. We will refine different numerical methods and try to solve the eigenvalue problem of the 
perturbation equation with high accuracy and efficiency. Comparing with other numerical methods, 
we find that the CFM and AIM methods can meet requirements of the high accuracy and efficiency in 
the computation. We will first review these two methods. Furthermore we will show that in different 
parameter ranges, the accuracy and efficiency also differ between these two refined methods. Without 
loss of generality, we will set M = 1 and fi = 0 in the following discussions. 
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3.1 Continued fraction method 


The CFM was proposed by Leaver m when he calculated the quasinormal modes (QNMs) of the 
Kerr black hole and is considered as the most accurate method to calculate the frequencies of pertur¬ 
bations m■ The core of this method is to cast the perturbation equation into a three-term recurrence 
relation, and from it we can get a continued fraction equation characterizing the perturbations. The 
CFM is thought to be able to give frequencies of perturbations with high numerical precision as there 
is no intermediate approximation compared to other numerical methods. See reviews HI ESI for more 
details. In this subsection, we try to get the three-term recurrence relation and the corresponding 
continued fraction equation. 

To calculate the frequency, we start from the physical boundary conditions which can be derived 
by studying the asymptotic behavior of Eq. Q 

( e~ t ^ u ~ 2 M') r *, r* —>• — oo(r —>• 2M), 

T ~ < —_ (8) 

I e *\A j2 -h 2? '*, r* —oo(r —>• oo), 


which means that there only exists the ingoing wave at the event horizon and the outgoing wave at 
the infinity. 

A solution to the radial equation Eq. (J6]) encoding the above boundary conditions can be written 
in the form as 


qj = 


r yi{2Mu-qQ) / r y2M(u,+V^ 2 -M 2 -#) iy /^Z^ {r _ 2M ) ^ „ n 

2 M ) \2MJ 6 ^ anX ’ 


(9) 


in which x = r and ao = 1. Substituting this expansion into the radial equation (0, we get a 
six-term recurrence relation. 


A)Oo + atocii = 0, 

7io 0 + Aiai + aia 2 = 0, 

A a o + 72ai + P 20-2 + 0:203 = 0, 

7300 + 5 3 ai + 7302 + /? 3 03 + 0304 = 0, 

04Oo + 7401 + <54d2 T 7403 + /?404 + 0:405 = 0, 

0nOn—4 T 7»iOn— 3 T Sn&n— 2 T r JnO’n— 1 T T Ot n a n -\-\ — 0, 


( 10 ) 
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where the recurrence coefficients are given by 


a n = ( Q 2 — 2) 2 [—4n 2 — 8(1 + iqQ — 2 ioj) n — 4(1 + 2 iqQ) + 16 iuS\ , 

P n = 2(Q 2 — 2) [2n 2 (5Q 2 — 6) + 4n (—2 + iqQ(5Q 2 — 6) + Q 2 (l — 14iu;) + 20iw)] (11) 

+2(Q 2 — 2) [—4 — 41(1 + 1) + Q 2 + 4 qQ{i — qQ)(Q 2 — 2)] 

+2(Q 2 - 2) [8(5gQ - 2 i)(Q 2 - 2)w - 64(Q 2 - 2)w 2 ] , 

7n = -8n 2 (5Q 4 -12Q 2 + 6) + 16/(/ + l)(Q 2 -l) 

-16 + 36Q 2 - 17Q 4 + 16iqQ(2 - 6 Q 2 + 3Q 4 ) + 4g 2 Q 2 (12 - 20Q 2 + 7Q 4 ) 

-16 [*(8 - 20Q 2 + 9Q 4 ) + 4gQ(6 - 11Q 2 + 4Q 4 )] u + 64(12 - 20Q 2 + 7Q 4 )cj 2 
-16m [2 i + qQ {6 - 12Q 2 + 5Q 4 ) + Q 4 (3i - 17w) - 24w + Q 2 (44w - 6z)] , 

6 n = 8n 2 (2 — 8Q 2 + 5Q 4 ) — 81(1 + 1)Q 2 (12) 

+16 - 100Q 2 + 79Q 4 - 16iqQ(2 - 10Q 2 + 7Q 4 ) - 4q 2 Q 2 (4 - 16Q 2 + 9Q 4 ) 

+16 [27(4 - 20Q 2 + 13Q 4 ) + qQ(8 - 32 Q 2 + 19Q 4 )] u - 64(4 - 16Q 2 + 9Q 4 )w 2 
+16m [2 i + qQ {2 - 8 Q 2 + 5Q 4 ) + Q 4 (7i - 19w) - 8 lu + 2Q 2 (16w - 5 i)] , 

rin = Q 2 [—4n 2 (5Q 2 - 4) + 4q 2 Q 2 (5Q 2 - 4)] 

+8nQ 2 [-8 - iqQ(5Q 2 - 4) + Q 2 (ll + 20 iu) - 16 iu\ (13) 

+Q 2 (-5i + 8 u) [12* - 32w + Q 2 (-19i + 40a;)] 

+8 qQ 3 [87 - 16w + Q 2 (-117 + 20a;)] , 

9 n = Q 4 [4n 2 + 87n(37 + qQ — 4 oj) + 35 — 4 q 2 Q 2 + 967a; — 64 cj 2 + 8qQ(—‘3i + 4a;)] . 


We can use the Gauss elimination to reduce the six-term recurrence relation to a five-term recur¬ 
rence relation, 


Po = A),7i = 7i, Pi = Pi 

73 = 73^3 = ^ 3,73 = 73 

/ 9 n ,-/ 

7 n 7 n / Vn —1 

7n-l 

/ On qi 

7 n 7 n “7 Pn—1 


, s' 2 = = 72,^2 = P 2 , 

, #3 = ft+J, = a n (« > 0 ), 

, <5n = i(« > 4 ), 

7n-l 

> ^ = ^-T^aLi(n>4). 
7 n —1 


( 14 ) 
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By repeating the Gauss elimination, a four-term recurrence relation can be derived, 



( 15 ) 


And at last we get a three-term recurrence relation 



(16) 


Then the frequencies of the perturbations are the solutions to the characteristic continued fraction 
equation 



(17) 


It is obvious that the six-term recurrence relation reduces to a four-term recurrence relation when 
Q = 0 or Q = \/2 due to the vanish of some coefficients in (1121) . When Q = 0, the GHS metric reduces 


to the Schwarzschild metric, and the six-term recurrence relation boils down to a four-term recurrence 
relation which later coincides with the three-term recurrence relation derived in ref. m after doing 
the Gauss elimination with 2 M = 1. 

3.2 Asymptotic iteration method 

The asymptotic iteration method (AIM) was first used to solve the eigenvalue problems of second 
order homogeneous linear differential equations [16]. It was then applied to find the frequencies of 
perturbations in Schwarzschild and Schwarzschild (anti-) de Sitter black holes |l7j . See review [18] 
and references therein. 

Let’s consider a second order homogeneous linear differential equations 


X" = ^o(x)x' + s 0 (x)x, 


(18) 


where Ao(x) and So(x) are smooth functions in some interval [a, b\. Differentiating it with respect to 
x, we get 
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where 


Ai(x) = A' 0 (x) + s 0 (:c) + (A 0 (.t)) 2 , si = Sg(x) + s 0 (x)A 0 (x). 
Using this step iteratively, we can get the (n + 2)th derivatives 

X (n+2) = K(x)x' + S n (x)x 


( 20 ) 


( 21 ) 


where 

A n (x) = A! n _x(x) + s n _i(x) + Ao(x)A n _i(x),s n (x) = 4 -iOe) + s 0 (x)A n _i(x). (22) 

For sufficiently large n, the asymptotic aspect of the method was introduced |18j . 

s n (x) s n - i(x) 


Art(^-) A n _i(x) 


(23) 


which is equivalent to imposing a termination to the number of iterations [23J. The perturbation 
frequencies can be derived from this “quantization condition”. However, the derivatives of A n (x) and 
s n (x ) in each iteration slow down the AIM considerably and also lead to precision problems. These 
drawbacks were overcomed in ref. HZ]. One can expand A n (x) and s n (x ) in Taylor series around a 
regular point £ at which the AIM is performed, 


-MO =^2 c h( x -tY i Sn(x) =^d' n (x-t)\ 

i=0 


(24) 


i=0 


Here c l n and d l n are the ztli Taylor coefficients of A n (£) and s n (£), respectively. Substituting these 
expansions into (1221) . we get a set of recursion relations for the coefficients, 


c n — (* + 1 )° l n\ + dn- 1 + C 0 C n-l > dn — (i + ^)d l +\ + ^ ^0 c ra-l- 


k=0 


k=0 


The quantization condition then can be expressed as 

ffi _ jO o n 
u ?i L n—1 u n—l L n 


(25) 


(26) 


which will give us the perturbation frequencies of a black hole. Both the accuracy and efficiency of 
AIM are greatly improved with this expansion [El- 

Now we apply this method to the GHS black hole. Taking a coordinate transformation x = 1— 
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and an abbreviation a 



the perturbation equation ([6]) turns into the standard form as Eq. m 


d 2 'S> 

_ = AoW _+ S|)W * 


(27) 


in which 


A 0 (x) 

s 0 (x) 


3x — 1 
x(l — x) ’ 


-1 


x 2 (1 — x) 2 

1 

x 2 (l — x) 2 


3 / 2 — a + ax 

4 \ 1 — a + ax 

f 2 Muj 


x 2 — 


1 — a + ax 


y 


1 — x 


— qQ — (1 — a + ax) 


1 

2 

4M 2 x 


1 — a + (a + l)x 
1 — a + ax 

1(1 + l)x 


(l-x) : 


rAi 


(1 — a + ax) 


(28) 

2 — a + ax 
1 — a + ax 



Now the infinity corresponds to x —> 1 and the horizon is at x 0. We can choose the regular point 
£ between 0 and 1. Substituting (1281) into the (1241) (1251) (1261) . we can get frequencies of perturbations 
around the GHS black hole. 

As we can see from (124[) . the efficiency and accuracy of the numerical result depend on the position 
of the expansion £. In our calculation, we find that when the charge of the black hole is not too large, 
the position of the expansion point has little influence and the AIM converges well. However, when 
the charge of the black hole is larger than M, the position of expansion will affect the efficiency and 
accuracy of the numerical computation apparently. 


3.3 The efficiency and accuracy comparisons between AIM and CFM 

In this subsection, we present the efficiency and accuracy comparisons between two numerical methods, 
the CFM and the AIM. The reason for us to concentrate on these two numerical methods in doing 
computation is that we have found that other numerical methods, such as the shooting method |241125] . 
the WKB method [261 1221 :2SJ and the finite difference method [2S] G2I 21] cannot give us good 
convergence and reliability in the computation when the scalar field is charged, although they can give 
consistent frequencies for the neutral scalar perturbations. In our numerical computations we have 
set M = 1 and (i = 0. 

For the fundamental modes of the charged scalar perturbation, it is easy to see that both methods, 
the AIM and the CFM, are easy to converge. For example as shown in Table 1, when l = 1,Q = 0.5 
and q = 1, only 20 iterative steps are needed for CFM to get the fundamental mode with relative error 
10~ 5 compared to the result obtained using 100 iterative steps, and 30 iterative steps for AIM with 
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relative error 10~ 4 compared to the results obtained by 100 iterative steps, respectively. However, for 
the overtones, we found that more iterative steps are needed compared to the fundamental modes in 
order to keep the accuracy. We adopt 100 iterative steps for both AIM and CFM in our numerical 
calculations. 


l = 1,Q = 0.5, q = 1 

CFM 

AIM (£ = 0.45) 

iterative steps 

20 

30 

fundamental n = 0 

(relative error with respect to 100 steps) 

0.491714-0.111583i 

(0.0005%,0.0028%) 

0.491758-0.111572i 

(-0.0085%,0.0128%) 

overtone n = 1 

(relative error with respect to 100 steps) 

0.479728-0.339701i 

(0.0194%,-0.0126%) 

0.477759-0.339279 

(4.2983%,0.1115%) 


Table 1: Fundamental mode n = 0 and the first overtone n = 1 when l = 1,Q = 0.5 and q = 1. The two numbers 
in brackets are relative errors of the real and imaginary part of the modes compared to the ones calculated with 100 
iteration steps, respectively. For the AIM, we do the expansion at the point £ = 0.45. 

When the charge Q of the black hole is small, results given by two methods, the CFM and AIM, 
agree quite well with each other. However, the situation changes when Q becomes large. We find that 
for the chosen large value of Q, the speed of convergence of AIM is faster than that of the CFM when 
the scalar perturbation field is weakly charged, but when the charge q of the scalar field increases, the 
speed of the convergence of the CFM becomes faster than that of the AIM. See Table 2 for concrete 
examples. 
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l = 1, n = 1 

in 

o 

II 

q = 2 

iterative steps 

80 

(relative error with respect to 100 steps) 

80 

(relative error with respect to 100 steps) 

Q = 0.1 

AIM 

(£ = 0.45) 

0.28448-0.31036i 
(-4.1 x 10 —6 , —3.0 x 10 -6 ) 

0.343659-0.321025i 
(-2.7 x 10 —6 , —1.4 x 10 -7 ) 

CFM 

0.28448-0.31036i 
(-2.1 x 10 -7 ,4.0 x 10~ 6 ) 

0.343658-0.321024i 
(1.3 x 10 —6 , 3.4 x 10“ 6 ) 

<5 = 1 

II > 
P £ 

Or 

0.54590-0.34803i 
(6.2 x 10 —8 , —4.1 x 10~ 7 ) 

1.198949-0.380392i 
(7.5 x 10 —6 ,1.6 x 10“ 5 ) 

CFM 

0.54590-0.34803i 
(-1.4 x 10 —6 ,3.2 x nr 7 ) 

1.198958-0.380398i 
(-1.7 x 10 —11 ,3.6 x 10" 11 ) 


Table 2: Convergence of CFM and AIM as parameters Q and q vary. When Q is small, Q = 0.1 for example, results 
given by both of these two methods have high accuracy and agree quite well. However, when Q becomes large, Q = 1 
for example, the convergence of the two methods depends on the value of q. For small q, q = 0.5 for example, the speed 
of the convergence of the AIM is faster than that of the CFM. However, as q increases, the speed of the convergence 
of the AIM slows down, while the speed of the convergence of the CFM increases and becomes faster than that of the 
AIM. 


We have made comparisons of these two methods in various situations, and the results are sum¬ 
marized in Table 3. The AIM and the CFM have little difference in the speed of convergence and 
accuracy for small Q. For large Q (but not larger than 1), the AIM converges faster than the CFM 
when q is small but when q is large we find that the CFM is a better way in computation. When Q 
is larger than 1, the CFM fails to converge and we can only rely on the AIM. 


The speed of 

convergence 

Small Q 

large Q 

AIM CFM 

AIM 

CFM 

Small q 

Little difference 

Faster 

Slower 

Large q 

Little difference 

Slower 

Faster 


Table 3: The speed of the convergence comparisons between the AIM and the CFM for calculating the overtone n = 1 
with various Q(Q < 1) and q. 

Moreover, it is known from (j24[) that the speed of convergence of AIM is related to the position of 
the expansion point. When Q is small, the position of the expansion £ lies in a large range to permit 
the convergence of the AIM. However as the increase of Q , the range of the allowed expansion point 
to accommodate convergence becomes narrower and the value of £ needs to be taken smaller. We 
choose £ = 0.45 when Q < 1 and £ = 0.43 when Q > 1 in our calculation. See the Table 4 for the 
comparison. For Q > 1, fine tuning of the expansion point is needed to control the convergence of the 
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computation by using the AIM. 


l = 2,Q= 1.2, q = 2 

AIM (£ = 0.45) 

AIM (£ = 0.43) 

iterative steps 

80 

80 

overtone n = 2 

(relative error with respect to 100 steps) 

1.695567-0.631905i 
(-4.7 x 10 -5 ,3.3 x 10” 4 ) 

1.695564-0.631910i 
(-3.9 x 10 —8 ,5.1 x 10” 7 ) 


Table 4: The dependence of the speed of the convergence in the AIM on different chosen positions of the expansion. 


4 Numerical results 

In this section we report the frequencies of the charged scalar perturbations in the stringy black holes 
with the change of the parameters, such as the charge of the black hole Q , the charge of the perturbing 
scalar field q and the angular momentum index l. We will analyze the numerical results to see their 
effects on the frequencies of the perturbation. 

4.1 The fundamental modes 

In this subsection, we study the frequencies of the fundamental modes of the charged scalar perturba¬ 
tions. In Table [5l we give the results for various Q and q when l = 0. We use both of the two highly 
precise methods to study the frequency domain of the perturbation, and find that they give consistent 
results when Q < 1. When Q > 1, for example Q = 1.2 in our table, we find that the CFM fails 
to converge so that we can only rely on the AIM in the computation. This challenges the argument 
that CFM is the most accurate method in calculating the frequencies of perturbations [19 1|. In our 
computations, we fix the number of iterative steps to be 100 in the two methods. 

We found that in all cases, the imaginary frequencies of the perturbations are negative, which 
indicate that there are no unstable modes of the charged scalar perturbation around the GHS black 
hole. When the black hole charge disappears, we recover the Schwarzschild black hole and the fun¬ 
damental modes we get in this limit reproduces the result of the Schwarzschild black hole under the 
scalar perturbation. Fixing the charge of the perturbation held, we find that the perturbation presents 
more oscillations but decays faster with the increase of the black hole charge Q. Compared with the 
real part, the imaginary part of the frequency changes slower when we vary Q or q. 

The objective pictures of the dependence of the real and imaginary parts of the frequencies on 
varying Q and q are shown in Figjl] From the left panel, we see that for the fixed Q, ujr increases 
linearly as q increases with a slope ~ ^ Q . We fitted the data and found that ujr ~ 0.079 + 0A93Qq. 
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l = 0 

0 = 0.1 

Q = 0.5 

Q = 1 

0 = 1.2 

g = 0 

0.11064-0.10493i 

0.11562-0.10599i 

0.13736-0.10961i 

0.15948-0.11120i 

q = 0.5 

0.13026-0.10804i 

0.21716-0.11732i 

0.34938-0.12550i 

0.42072-0.12873i 

g=l 

0.15019-0.11070i 

0.32485-0.12258i 

0.57549-0.12942i 

0.69534-0.13332i 

q = 1.5 

0.17044-0.11297i 

0.43710-0.12495i 

0.81000-0.12991i 

0.97805-0.13401i 

9 = 2 

0.19098-0.11491i 

0.55260-0.12596i 

1.04946-0.12948i 

1.26577-0.13343i 

9 = 2.5 

0.21180-0.11658i 

0.67038-0.12634i 

1.29192-0.12885i 

1.55668-0.13252i 

9 = 3 

0.23288-0.11800i 

0.78977-0.12643i 

1.53630-0.12827i 

1.84974-0.13159i 

9 = 4 

0.27578-0.12025i 

1.03171-0.12630i 

2.02853-0.12735i 

2.43985-0.13001i 

9 = 5 

0.31955-0.12188i 

1.27626-0.12607i 

2.52344-0.12673i 

3.03318-0.12886i 

9 = 6 

0.36409-0.12306 

1.52236-0.12587i 

3.01986-0.12632i 

3.62839-0.12803i 

9 = 7 

0.40928-0.12392i 

1.76944-0.125711 

3.51721-0.12603i 

4.22478-0.12743i 

9 = 8 

0.45505-0.12453i 

2.01718-0.12559i 

4.01517-0.12582i 

4.82197-0.12698i 

q= 10 

0.54797-0.12528i 

2.51392-0.125411 

5.01226-0.12555i 

6.01789-0.12637i 

9 = 12 

0.64236-0.12564i 

3.01169-0.12530i 

6.01027-0.12539i 

7.21530-0.125811 


Table 5: Fundamental modes of the GHS BH when l = 0. The CFM and the AIM give the consistent result when 
Q < 1. However, when Q > 1, q = 1.2 in the table for example, the CFM fails to converge and at this time we can only 
rely on the AIM. The iteration steps in both methods are taken to be 100. The expansion point in the AIM is taken to 
be £ = 0.43. For Schwarzschild BH ( Q = 0), the fundamental mode reads uj = 0.11047 — 0.10487i. 


The behavior of ujq is more complicated. From the right panel, we observe that when Q is small, 
Q = 0.1 for example, \ujj\ is a monotonically increasing function of q. However, when Q becomes 
large, |w/| is no longer a monotonic function of q. With the increase of q, |w/| first increases and 
then after reaching a maximum value it begins to decrease. With the increase of Q, we find that the 
imaginary part of the perturbation becomes more negative when the perturbation is weakly charged, 
which indicates that in this case the black hole is more stable. When q becomes big enough, w/ 
converges to a constant no matter what values of Q one chooses. 

We plot the effective potential in Fig|2]to give some intuitive understandings of the properties of 
the fundamental modes discussed above. We find that V e ff —» — fypr when r -A 2 M, while V e ff — > 0 
when r —» oo. When q = 1, the effective potential has a barrier for small Q. The barrier becomes 
lower and even disappears as Q increases. This implies that the perturbing wave can fall into the 
black hole more easily when Q becomes larger, which explains the faster decay of the weakly charged 
scalar perturbation with the increase of Q as shown in the right panel of Fig(T] On the other hand, 
near the horizon, the potential is more negative for larger Q , which tells us that the perturbation can 
have more momentum to fall into the black hole and explains the reason that the perturbing scalar 
can have faster oscillation in the decay when Q increases. This supports the observation in the right 


14 
























q 


q 


Figure 1: The fundamental modes of charged scalar perturbation of GHS BH when 1 = 0. 


panel of FigJT] that we have bigger real part of the perturbation frequency for bigger Q. 

When q becomes large enough, q = 10 for instance or even bigger, the potential barrier disappears 
for all Q. The perturbation wave can be absorbed by the black hole without any obstacles. This gives 
the same decay speed of the perturbation for all values of Q as shown in the right panel of FigQ] On 
the other hand, at the black hole horizon, the differences in the potential values there caused by the 
black hole charge Q become bigger compared to the case with weakly charged scalar perturbation. 
This explains that the difference in the momentum for the perturbing wave to fall into the hole is 
enlarged when the scalar field is heavily charged, which explains the big difference in the real part 
of the frequency with the change of the black hole charge Q when the scalar perturbation is heavily 
charged as shown in the left panel of Fig{T] 



r r 


Figure 2: The effective potential for different Qs when 1 = 0 (left for q = 1 , right for q = 10). 


We also perform the same calculations for different angular indexes l = 1 and l = 2, respectively, to 
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see the effect of the angular momentum. The results are listed in Tablc|6]and0 The overall behaviors 
of the fundamental modes for l = 1 and l = 2 are similar to that of l = 0. Also, we can see that the 
imaginary part of the fundamental modes changes slowly while the real part changes significantly as 
we vary Q and q. 


l = 1 

Q = 0.1 

Q = 0.5 

Q= 1 

Q = 1.2 

g = 0 

0.29343-0.09771i 

0.30622-0.09899i 

0.36304-0.103611 

0.42191-0.10624i 

q= 0.5 

0.31072-0.09940i 

0.39622-0.10623i 

0.55527-0.11470i 

0.66341-0.11785i 

9=1 

0.32829-0.lOlOOi 

0.49171-0.11158i 

0.76093-0.12112i 

0.91814-0.12457i 

q = 1.5 

0.34613-0.10250i 

0.59173-0.11549i 

0.97615-0.12478i 

1.18197-0.12850i 

9 = 2 

0.36422-0.103911 

0.69548-0.11833i 

1.19842-0.12680i 

1.45242-0.13075i 

9 = 2.5 

0.38256-0.10524i 

0.80232-0.12038i 

1.42597-0.12785i 

1.72788-0.13196i 

9 = 3 

0.40113-0.10649i 

0.91173-0.12185i 

1.65756-0.12832i 

2.00725-0.13252i 

9 = 4 

0.43895-0.10877i 

1.13661-0.12365i 

2.12944-0.12842i 

2.57463-0.13257i 

9 = 5 

0.47761-0.11077i 

1.36750-0.12457i 

2.60924-0.12808i 

3.15021-0.13200i 

9 = 6 

0.51704-0.11254i 

1.60265-0.12503i 

3.09417-0.12766i 

3.73137-0.13124i 

9 = 7 

0.55718-0.11409i 

1.84088-0.12525i 

3.58258-0.12726i 

4.31647-0.13048i 

9 = 8 

0.59799-0.11544i 

2.08137-0.12535i 

4.07342-0.12692i 

4.90445-0.12978i 

g= 10 

0.68138-0.11767i 

2.56706-0.12538i 

5.05994-0.12641i 

6.08632-0.12865i 

g = 12 

0.76686-0.11938i 

3.05688-0.12535i 

6.05055-0.12606i 

7.21341-0.12762i 


Table 6: Fundamental modes of the GHS black hole with the angular index 1 = 1. The CFM and the AIM give the 
same results when Q < 1. However, when Q > 1 ,q = 1.2 in the table for example, the CFM fails to converge and so 
that we can only rely on the AIM. The iteration steps in the two methods are taken to be 100. The expansion point in 
the AIM is taken to be £ = 0.43. For the limiting case (Q = 0), we reproduce the result of the Schwarzschild black hole 
with the fundamental mode uj = 0.29293 — 0.09766i. 

In Figj3] we plot the fundamental modes for different angular index l (l = 0,1, 2). From the figure, 
we can see that when Q is fixed, the real part of the frequency keeps linearly increasing with the 
increase of q, which is independent of the angular index l. For l = 1, ujr ~ 0.252 + 0.481 Qq. For 
l = 2, Wfl ~ 0.455 + 0A71Qq. Moreover, for fixed Q, ujrs for different Is approach to each other as q 
becomes very large. This is even clearer when the black hole charge Q is large. As q increases, the 
imaginary part uiq flattens and approaches to constant values regardless of the chosen Q and l. 

We plot the effective potential for different Is when Q = 0.5 in Figj4]to understand the curves in 
Fig|3] more intuitively. When q is small, the potential barrier increases as l increases. This means 
that the perturbing wave with higher l is more difficult to be absorbed into the black hole and decay 
so that its imaginary frequency is bigger. At the horizon, the potential drops faster when l is larger, 
thus for higher l, the perturbation falls into the black hole usually with bigger momentum so that it 
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can oscillate faster. At large q, the potential barriers disappear for all values of l. Furthermore, they 
basically coincide with each other for large q (q = 10 for example in Fig|3J. So both the real parts and 
the imaginary parts of the perturbation modes approach to each other at large q. as shown in Figj3] 



Figure 3: The fundamental modes of GHS BH for different Is and Qs (The first panel for real parts, the remaining 
panels for imaginary parts). 



4.2 Overtones 

Now we turn to present the results of the overtones of the charged scalar perturbation around the GHS 
black hole. We only list the fundamental modes and the first two overtones when l = 2 in Tabled for 
concision. 
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l = 2 

overtones 

Q = o.i 

Q = 0.5 

<3 = 1 

Q= 1.2 

O 

II 

n = 0 

n = 1 

n = 2 

0.48445-0.09681i 

0.46470-0.29575i 

0.43145-0.50876i 

0.50541-0.09812i 
0.48660-0.29942i 

0.45489-0.51401i 

0.59878-0.10291i 

0.58380-0.31271i 

0.55839-0.53288i 

0.69577-0.10572i 

0.68476-0.32025i 

0.66619-0.54286i 

q = 2 

n = 0 

n = 1 

n = 2 

0.55330-0.10086i 

0.53702-0.30696i 

0.50965-0.52437i 

0.87667-0.11296i 

0.86989-0.34051i 

0.85816-0.57223i 

1.39915-0.12269i 

1.39785-0.36838i 

1.39555-0.61487i 

1.69259-0.12634i 

1.69356-0.37907i 

1.69556-0.63191i 

q = 4 

n = 0 

n = 1 

n = 2 

0.62484-0.10437i 

0.61152-0.31664i 

0.58903-0.53793i 

1.29143-0.11985i 

1.28926-0.35996i 

1.28531-0.60114i 

2.28913-0.12753i 

2.28956-0.38258i 

2.29043-0.63758i 

2.77507-0.13167i 

2.77690-0.39496i 

2.78057-0.65810i 

5 = 6 

n = 0 

n = 1 

n = 2 

0.69884-0.10739i 

0.68798-0.32498i 

0.66956-0.54970i 

1.73335-0.12292i 

1.73272-0.36885i 

1.73153-0.61501i 

3.22242-0.12824i 

3.22267-0.38472i 

3.22316-0.64124i 

3.89915-0.13246i 

5 = 8 

n = 0 

n = 1 

n = 2 

0.77505-0.10998i 

0.76626-0.33215i 

0.75120-0.55992i 

2.19246-0.12427i 

2.19230-0.37283i 

2.19200-0.62143i 

4.17893- 0.12791i 

4.17894- 0.38374i 

4.17898-0.63962i 


O 

T—1 

II 

OH 

n = 0 

n = 1 

n = 2 

0.85328-0.112201 

0.84618-0.33831i 

0.83392-0.56879i 

2.66269-0.12486i 

2.66267-0.37459i 

2.66262-0.62432i 

5.14884-0.12741i 

5.14876-0.38224i 

5.14859-0.63710i 


5 = 12 

n = 0 

n = 1 

n = 2 

0.93333-0.11409i 

0.92763-0.34359i 

0.91766-0.57647i 

3.14034- 0.12511i 

3.14034- 0.37534i 

3.14036-0.62557i 

6.12702-0.12696i 

6.12691-0.38089i 

6.12669-0.63482i 



Table 7: The frequencies of the charged scalar perturbations around the GHS black hole when l = 2. The CFM and 
the AIM give the same result when Q < 1. However, the CFM fails when Q > 1. The frequencies for Q = 1.2 are 
calculated by the AIM. We adjust the expansion point to keep high precision and efficiency. The iteration step is 100. 
When Q = 1.2 and q > 6, the precision of the results is lower than 10~ 5 and the results are not shown in this table. For 
the Schwarzschild black hole (Q = 0), the fundamental mode reads oj = 0.48364 — 0.09675i. 


It can be seen from Tabled that when Q, q and l are fixed, the real parts of the overtones are 
nearly invariant while the imaginary parts change more significantly. To see more clearly, we plot the 
overtones for l = 1 and / = 2 in Fig. 5. We can see that the real part ujr is nearly independent of n, 
while the imaginary part uj depends on n linearly. For example, when l = 1, wj can be approximated 
well by the form ujj ~ 0.1262 + 0.2562n when q = 6 and ujj ~ 0.1261 + 0.2520n when q = 12. As q 
increases, the slope tends to 0.252. Meanwhile, the modes for different l also become closer to each 
other as q increases. This can be understood since the effective potential for different l approaches to 
each other when q is large. 
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Figure 5: Overtones of perturbations around GHS black hole when Q = 1. The dashed left red line corresponds to 
the overtones when 1 = 1 and the right green line corresponds to the overtones when 1 = 2. 


4.3 The superradiance and its stability 

In this subsection we will show that the superradiance with the amplification of the incident wave can 
not appear for the GHS black hole. We consider the classical scattering problem for a charged scalar 
field in the GHS black hole background. 

The asymptotic behavior of ([6]) when w 2 > /x 2 can be derived straightforwardly, 


\J/ rv 


iJbfl-n 2 




„• r 


r* —» —oo(r —2M), 
r* —>• oo(r —>• oo). 


(29) 


This boundary condition corresponds to an incident wave of unit amplitude from the infinity and a 
reflected wave of amplitude B back to the infinity and a transmitted wave of amplitude T towards the 
event horizon. Since the effective potential is real, T* is also a solution of the radial equation Eq. (J6|) 
and independent of 4'. Thus, the Wronskian W = — 4'*^-\k is independent of r*. Calculating 

the Wronskian at the black hole horizon and at the infinity respectively, and equaling the two values, 
we get 


1 - 


B\ 2 


oj — 


qQ 

2 M 




\ T\‘ 


(30) 
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We see that if n < lo < the amplitude of the reflected wave is larger than the one of the incident 
wave |E>| 2 > 1. This phenomenon is known as superradiance. Thus, we get the superradiance condition 




< oj < 


qQ_ 

2M 


(31) 


This condition was also derived in |321 [33] . 

Multiplying the complex conjugated field T* on both sides of (JUJ) and doing partial integration, we 
get 

/ oo n oo 

(w - q^f \^{r*)\ 2 dr* = / (H|T(r*)| 2 + |T'(r*)| 2 ) dr*, (32) 

-oo J — OO 

in which <f> = y. Note that the right hand side of the above equation is real and positive since V is 
positive. Taking imaginary part of both sides, we get 


1 ro o 

(a 2 + 6 2 ) 4 cos ( - arctan -J + ujr - q$ h + J 2ujj (ojr - q$) |T(r*)| 2 dr* = 0. (33) 

Here we have taken notations u = ojr + iuj, a = uj 2 r — uj — b = 2 urujj and Since 

arctan ^ G ( — f j §) an< ^ d> is a monotonic decreasing function, we get ojj < 0 if Ur > q<&h- Thus, if 
the instability with ujj > 0 occurs, we must have ujr < q&h- This tells us that the instability can only 
take place when 

H <u R < q$ h , (34) 

which is the superradiance condition Eq. (13111 . Thus if the GHS black hole experiences instability, this 
instability must be superradiant instability. 

But from the above computations, we always find that the perturbation modes with lor are beyond 
the superradiant condition Eq. (ED, no matter what parameter ranges we choose. According to 
the analytical argument, these modes should be stable with uij < 0, which are consistent with our 
numerical results above. On the other hand, from the effective potential, we do not see the potential 
well outside the black hole to accumulate the energy. The necessary condition for the existence of 
the superradiant instability does not hold. Thus the superradiant instability cannot exist in the GHS 
background. 
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5 Summary and discussion 


In this paper, we have discussed the stability of the GHS black hole under charged scalar perturbations. 
For the charged scalar perturbations, we have found that not all available numerical methods can 
efficiently compute the accurate frequencies of the perturbations. This is different from that of the 
neutral scalar perturbations. We have discovered that two numerical methods, the CFM and the AIM, 
can still keep high accuracy and efficiency in the computations of the frequency of the charged scalar 
perturbations. We have experienced that the speed of convergence of the AIM depends on the position 
of the expansion point which needs to be chosen suitably. Besides, we have observed that the speed of 
convergence of CFM is close to that of the AIM when the black hole charge Q is small. However, when 
Q becomes large but still smaller than the mass of the black hole M, the AIM has faster convergence 
than the CFM when the scalar field is weakly charged; but the result is opposite when the charge of 
the scalar field is big. When the black hole charge Q exceeds the black hole mass M, the CFM is found 
invalid to give reliable results, while the AIM can still work. However, the convergence of the AIM 
becomes bad as Q approaches to the extremal value y/2. The results have been summarized in Table 
3. The comparisons of the efficiency and accuracy among different numerical methods are important, 
because the accuracy and convergence of the numerical computations are key requirements to grasp 
the properties in the perturbations around black holes. 

The influences on the frequencies of the perturbations brought by parameters describing the back¬ 
ground and the perturbation field have been illustrated. The intuitive reasons behind these phenomena 
have also been discussed. We have observed that the GHS black hole spacetime, which can reduce to 
the Schwarzschild black hole when the black hole charge goes to zero, is stable against the charged 
scalar perturbation. This result is different from what have been disclosed in the AdS black holes and 
the dS black holes with vanishing angular momentum, where the backgrounds experience instability 
under the charged scalar perturbations, while keep stable against neutral scalar perturbations. 

In |7| it was concluded that the new instability in the dS black hole against charged scalar pertur¬ 
bations with vanishing angular momentum is caused by the superradiance. This was further confirmed 
in [8]. For the GHS black hole spacetime, we have shown that the superradiance does not happen. 
Recently the superradiant instability of the GHS background was discussed in |34l 135 ], But they 
put the black hole in the artificial cavity. Although in this way they claimed that they can devise a 
black hole bomb, the mechanism with the artificial mirror is not convincing. One needs to introduce 
a natural wall, for example the massive fields [Ml EH, to trigger the superradiant instability. It was 
proved in |32| that the reflecting mirror made by the mass term of the incident field cannot trigger the 
superradiant instability and create the GHS black hole bomb. Our result has further argued that the 
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GHS black hole background is always stable against the charged scalar perturbation. Considering that 
the incident charged scalar wave cannot be amplified due to the superradiance around the GHS hole 
and the property of the effective potential, we find that the superradiant instability cannot happen in 
the GHS black hole background. 
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